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ABSTRACT 


The predictions of the two-phase accretion disc- 
corona models for active galactic nuclei are compared 
with observations. We discuss the possibility to use 
X-ray spectral slopes, equivalent widths of the iron 
line, and the observed flux-spectral index correlation 
as diagnostics of the X- 7 -ray source compactness and 
geometry as well as of the cold disc temperature. As 
an example of the application of the modelling tools, 
we use XSPEC to fit the broad-band data of Seyfert 
I galaxy, IC4329A, with a theoretical spectrum from 
a hemisphere-corona. 
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1. INTRODUCTION 


Observations of Seyfert 1 galaxies by Ginga and ear¬ 
lier by HE AO 1 showed that their X-ray spectra con¬ 
sist of at least two distinct components: an intrin¬ 
sic power-law with energy index, 02-18 =0.95±0.15, 
2 - 18 keV range and a Compton reflection 


bump (jM ushotzkv. Done, fc Po unds 1993 , Nandra 
fc Pounds 1994 [weaver et al. 199^). Becent ASCA 


obse: 


vations also show these features (|Nandra et al. 


19961) . OSSE obser vations show a much steeper spec- 


trum with a w 1.5 (Johnson et al. 1996]), an d COMP- 
TEL has no t detected Seyfert galaxies at all ( Maisack 
et al. 19^) . Averaged broad-band spectrum tor a 
sample of Seyfert 1 galaxies which was obtained us¬ 
ing non-simultaneous Ginga and OSSE data shows 
a high-energy cutoff of the intrinsic power-law at 
Ec =600^300 keV ( Zdziarski et al. 1995|). A sim ilar 


cutoff energy is obtained by pondek et al. (1996) 


us¬ 


ing the data from EXOSAT and USSE. HUSAT, Ginga 
and OSSE data of the Seyfert 1 galaxy, IC4 329A. also 


require a q utoff in the intrinsic spectrum (Madejski 
et al. 1995 ). The brightest Seyfert galaxy m 7 -rays, 
NCC 4151, shows a cutoff at ~ 150 keV, but is prob- 
^Iv a freak object re quiring special consideration s 
(iPoutanen et al 1996| , but see |Zdziarski et al. 19^ ). 

Non-thermal pair models intensively studied in the 


mid 80s (pvensson 1986, for review see Svensson 


1994 ) predicted the spectral index of the intrinsic 
: 0.9 — 1.0 and a prominent pair 
Absence of the annihilation line 


component Oint 
annihilation line, 
in the observed spectra rules out pure non-thermal 
models. Attention is now focused on thermal models. 
A definite answer on the question about the relative 
fraction of the thermal and non-thermal plasma can¬ 
not be obtained without high quality data above 200 
keV. 



Eigure 1: An example of the deconvolution of the emer¬ 
gent spectrum from a disc-corona system into differ¬ 
ent scattering orders. Solid curve represents the total 
emergent spectrum. Dotted curves give the contribu¬ 
tion of the unscattered radiation which consists of three 
components: black-body radiation with a maximum at 
€ = hvImeC? a Compton reflection compo¬ 

nent with a maximum at e ^0.1, and a broad annihila¬ 
tion line at e ~. 8 . Dashed curve gives the contribution of 
the single scattered radiation (notice that the broad bump 
with a maximum at e ^0.2 is the first scattering order of 
the Compton reflection component). All other scattering 
orders are represented by dash-dotted curves. 

There is a consensus that the X/ 7 -ray spectrum of 
Seyferts is produced by Comptonization in hot plas¬ 
mas of soft radiation from the UV source. The exact 
geometry of both phases is unknown. Presently, the 
most commonly use d scenario is the two-phas e disc- 


corona model (e.g., [Haardt fc Maraschi 1991, 1993) 


in which a hot X-ray emitting corona is located above 
the cold UV-emitting disc of the canonical black hole 



















































model for AGNs. The power-law X-ray spectrum 
is generated by thermal Comptonization of the soft 
UV-radiation. About half of the X-rays enters and 
is reprocessed by the cold disc, emerging mostly as 
black body disc radiation in the UV (some fraction 
is reflec ted producing a Compton r eflection compo¬ 
nent). Haardt & Maraschi (1991) emphasized the 


coupling between the two phases due to the repro¬ 
cessing, as the soft disc photons influence the cooling 
of the corona. They showed that nearly all power 
must be dissipated in the corona in order to have 
aint ~ 0.9-1. A consequence of this is that the soft 
disc luminosity, Lg, is of the same order as the hard 
X-ray luminosity, Lh. 


ported in|Stern et al. 1995b|. Two methods were used. 
The first method is based on the Non-Linear Monte- 
Carlo method bv S tern (see detailed description in 
Stern et al. 1995a). The second method is the it- 
erative scattering method (ISM) where the radiative 
transfer is exactly solved for each sca ttering order 
separately (Poutanen & Svensson 1996). The results 
of both codes are m excellent agreement. In this pa¬ 
per we report the results using the second method. 
A typical spectrum emerging from the disc-corona 
system and its different components is shown in Fig¬ 
ure ||. 



2. SPECTRA FROM THE DISC-CORONA 
SYSTEM 


The ISM code is a ID code but it can also treat quasi- 
2D radiative transfer in cylinders and hemispheres. 
The full Compton scattering matrix is used allowing 
us to treat polarized radiative transfer in thermal rel¬ 
ativistic plasmas. Fully angular dependent, polarized 
Compto n reflection is implemented using a Gr een’s 
matrix (Poutanen, Nagendra, & Svensson 19961) . For 
a given geometry, energy balance gives a unique 
combination of electron temperature 0 = kTe/meC^ 
and Thomson optical depth, tt = UeCTh. Solv¬ 
ing the pair balance for a given (0, ttI provides 
a unique diss ipation compactness, Idiss (phisellini & 
Haardt 1994). Here, the local dissipation compact- 
ness, Idiss = {Ldiss/h){aT/meC^), characterizes the 
dissipation with Ldiss being the power providing uni¬ 
form heating in a cubic volume of size h in the case 
of a slab of height h, or in the whole volume in the 
case of an active region (cylinder or hemisphere) of 
height h. The relation between 0, tt. and for 


different geometries is g iven in Stern et al. (1995b) 
(see also Svensson 199'*' 


Figure 2: Upper panel: speetra emerging in different di¬ 
rections from a slab-corona. Solid, dotted, and dashed 
curves represent the flux at viewing cosine angle pL=0.9, 
0.5 and 0.1, respectively. The temperature of the black- 
body radiation from the accretion disc is fixed at T\,\y=50 
e V. The electron temperature of the corona increases from 
the top of the figure to the bottom: 0 = kTe/mc^ = 
0.16, 0.2, 0.25, 0.32, O.f, 0.5, 0.63, 0.79, 1.0. Notice 
the increase of the iron line eguivalent width with increas¬ 
ing temperature due to anisotropic scattering effects (see, 
e.g., Haardt & Maraschi 1993, Stem et al. 1995b). The 
iron line profile has a triangular shape due to the energy 
resolution employed here. Lower panel: spectra from a 
cylinder-corona with a height-to-radius ratio 2. 


Observations show that Lg may be several times 
larger than Lh, in contradiction to the prediction of 
the uniform two-phase disc-corona mo del. This led 


Haardt, Maraschi, & Ghisellini (1994) to propose a 
patchy disc-corona model, where the corona consists 
of several localized active regions on the surface of the 
disc. Internal disc dissipation results in UV-radiation 
that leaves the disc without entering the active re¬ 
gions, thus enhancing the observed Lg relative to Lh. 


The first exact radiative transfer/Comptonization 
calculations in a disc-corona system accounting for 
energy and pair balance as well as reprocessing by 
the cold disc (including angular anisotropy) were re- 


Figure g shows spectra from slab-corona and 
cylinder-corona at different 0. As the temperature 
increases the first order scattered photons are pre¬ 
dominantly scattered back into the disc, causing a 
deficiency of the first order photons in the observed 
spectrum. This anisotropic eff ect gives rise to a n 
anisotropy break in the spectra (Stern et al. 1995b|). 


3. DIAGNOSTICS USING SPECTRAL SLOPES 
AND IRON LINE EQUIVALENT WIDTHS 


In order to compare the predictions of the two-phase 
corona model with observations, we determine the 
least square overall spectral slope, a 2 -i 8 : for theo¬ 
retical model spectra and display them in Figure ^ 
as a function of the dissipation compactness, ldiss, for 
different geometries. The right panels show the ob¬ 
served distributions of the spectr al indices for Ginaa 


spect ra of 27 Seyfert galaxies (Nandra fc Pound; 


1994) and ior ASCA spectra f2 - iU iteV range) ot 


15 beyfert Is (pMandra et al. 1996|). One sees that the 
observations are more consistent with active surface 
regions (such as hemispheres or cylinders) than with 
slabs. Active regions produce spectra covering the 
observed ranges of spectral indices and cutoff energies 
for the observed range of compactnesses. 
























































For small temperature of the black-body, Tbb, 
anisotropic effects significantly change the spectral 
index at compactnesses smaller than ~ 100. For 
2bb=50 eV, the anisotropy break in the 2-18 keV 
range appears all the way up to Idiss ~ 1000. For 
such a high temperature of the soft radiation, the 
first scattering order extends up to the Ginga range 
already for ~100 keV. 


When the anisotropy break shifts to the spectral re¬ 
gion around the fluorescent iron line at 6.4 keV, the 
continuum flux in the direction of the cold disc ex¬ 
ceeds the flux in the observer’s direction, causing an 
increase of the equivalent width of the K-line in com¬ 
parison with the isotropic case. This effect is demon¬ 
strated in Figure |4 where the iron line equivalent 
width is plotted against the electron temperature for 
different geometries. We also show the dist ribution 
of the equi valent widths observed by ASCA ( Nandra 
et d. 19961) . 
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Figure 3: Overall power-law index, 02-18, least- 

square fitted to the model spectra in the 2 - 18 keV 
energy range vs. dissipation compactness, l^iss = 
(idiss/h)(crT/4neC^)- Rectangles, hemispheres, and 
squares represent results for the slab-corona, hemisphere- 
corona, and cylinder-corona (with height-to-radius ratio 
2), respectively. Filled symbols correspond to the spectra 
m almost face-on direction (fi—0.9), open symbols rep¬ 
resent results for inclination 6CP (pL=0.5). Results for 
Thh=5 eV and T\^\,=50 eV are presented in upper and 
lower panels, respectively. Lower compactness limits rep¬ 
resent Seyfert galaxies that have known estimates of their 
X-ray variability. At, and thus lower limits of their com¬ 
pactnesses, lobs = iLo\,s/cAt){aT/rrieC^), see Done & 
Fabian 1989. The right panels shows the observed dis¬ 
tributions of power-law indices of Seyfert galaxies. Solid 
histogram: Ginga data from Nandra & Pounds (1994). 
Dotted histogram: ASGA data from Nandra et al. (1996, 
Table 2). 



Figure 4' Equivalent width of the iron K-line at 6.4 keV 
vs. electron temperature, 0. Same notations as in Fig¬ 
ure 3. Notice that for low temperature, the equivalent 
width is almost constant. R increases when the tempera¬ 
ture increases due to the anisotropy effects (see Fig. 2). 
For higher disc temperature, Tbb, the anisotropy effects 
become important at smaller electron temperature. Right 
panels show the distribution of the equivalent widths of 
Seyfert Is observed by ASGA (Nandra et al. 1996, Ta¬ 
ble 6). 



Figure 5: Gorrelation between spectral index, 02 - 13 , o,i^d 
the continuum flux in 2 - 10 keV energy band corrected for 
absorption (in units of 10~^^ erg s“^ cm“^^ observed by 
Ginga in NGG 4151. Filled spheres with the error bars 
show the data from Yaqoob & Warwick (1991) and open 
spheres are the data from Yaqoob et al. (1993). Hemi¬ 
spheres and solid curve show the correlation between spec¬ 
tral index and dissipation compactness (arbitrary com¬ 
pactness units) for the hemisphere-corona model. Similar 
correlation is expected also in non-thermal pair models 
(see Yaqoob 1992). 




























































4. FLUX - SPECTRAL INDEX CORRELATION 


Some Seyfert galaxies show correla ted variations in 


flux a nd spectral index f see. e.g.. Matsuoka et al . 
IQilol tTreves et al. 1990L |Yaqoob & Warwick iyyi|. 
YaqooB^et al. 19^). This correlation has been ex¬ 


plained in the context of non-thermal pair model (e.g. 
Yaqoob 1992|). Our thermal pair model shows strong 
correlation between compactness parameter and the 
spectral index (see Fig. Assuming that the size 
of the X /7 source is con^ant there is a direct corre¬ 
spondence between /diss and the observed flux (here 
we neglect the eff^t of spectral changes with com¬ 
pactness). Figure g shows the flux - spectral index 
correlation observed in NGC 4151 and the compact¬ 
ness - spectral index correlation for a hemisphere- 
corona (^diss is in arbitrary units and is rescaled to 
fit observations). It can be seen in Figure^, that for 
Tbb=5 eV spectral index changes dramatically from 
1.0 to 0.5 when Zdiss changes by an order of magnitude 
from 100 to 10 . 


to set constraints on the viewing angle within the 
framework of the two-phase pair corona model, once 
high quality spectra become available. (2) The spec¬ 
tra also depend on the geometry of the coronal re¬ 
gions, so observed high quality spectra can be used 
as diagnostics of the geometry. Presently, it seems 
that active regions are favored over homogeneous slab 
coronae. (3) Observations of flux-index correlations 
can provide constraints on the source compactness. 
(4) The equivalent width of the iron line can be 
used as a diagnostic of the temperature of the repro¬ 
cessed (black-body) radiation which then constrains 
the size. 
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5. SPECTRAL FITTING 


We implemented the ISM code into XSPEC, and used 
it to fit the broad band spectrum of the second bright¬ 
est Seyfert 1 galaxy, IC4329A. As an example, we fit¬ 
ted data with theoretical spectra from a hemisphere- 
corona. The spectrum is absorbed by a combination 
of neutral and ionized absorbers. It is presently difR- 
cult to constrain the temperature of the black-body 
radiation, Tbb, and hence we fixed it at 10 eV. The 
best fit model parameters are Te= 100 l'l 3 Q keV, and 
TT= 1 . 3 to; 2 ' The quality of the data does not allow 
us to constrain the geometry of the emitting region. 
The data and the model spectrum are shown in Fig¬ 
ure n 



E ( keV ) 

Figure 6: The broad-band spectrum of Seyfert 1 galaxy, 
IC 4329 A. ROSAT, Ginga, and OSSE observations are 
described in Madejski et al. (1995). The upper limits 
are 2a. Solid curve shows the best fit hemisphere-corona 
model spectrum. 


6 . CONCLUSIONS 


Several diagnostics from physical modelling of ob¬ 
served spectra are possible. (1) As anisotropic ef¬ 
fects are very important, i.e. as spectral shapes de¬ 
pend strongly on viewing angle, it will be possible 
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